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Abstract 

We have searched for the leptonic and radiative leptonic B decays, — > ii^ffx, B~^ — > e'^Vel 
and S+ jjL^v^^. Using a 140 fb~^ data sample collected with the Belle detector at the KEKB 
asymmetric e'^e" collider, we find no evidence for signals in any mode and set the following 
preliminary upper limits at 90% confidence level: B{B'^ — ^ l^^^n) < 2.0 x 10""^, B{B'^ e'^u^j) < 
2.2 X 10-5 and B{B+ m"^ 1/^,7) < 2.3 x 10"^ 

PACS numbers: 13.25.Hw, 11.30.Er, 12.15.Hh 
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I. INTRODUCTION 



In the Standard Model (SM), the fully-leptonic decays i^i^e, where i represents an 

electron, muon, or r lepton, are allowed through annihilation into a virtual W boson. The 
branching fraction is given by: 
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where Gi;' is the Fermi coupling constant, itlb and are the masses of the B meson and 
lepton, tb is the B meson lifetime, Vub is an element of the Cabibbo-Kobayashi-Maskawa 
quark mixing matrix 0, and is the decay constant that parameterizes the overlap of 
the quark wave functions within the meson. A measurement of this branching fraction, 
combined with the value of \Vub\ obtained from other decay modes such as i? ^ 7r£z7£, allows 
us to determine /b, which is needed to extract \Vtd\ from measurements of B^B'^ mixing. 
The theoretical expectation for the branching fraction of B^ —>■ t^Vt lies in the range 
(1 — 10) X 10~^. Since the branching fraction is proportional to the square of the mass of the 
charged lepton, the branching fractions for B^ f^^^^fi and B^ —>■ e^Vf. are suppressed by 
factors of 225 and 10^, respectively {^^helicity suppression"). Although branching fractions 
at these levels are beyond the reach of the Belle experiment with the current data sample, 
an observation of either of the latter decay modes would be a clear indication of physics 
beyond the SM. For example, the B^ £"*"z/^ decay rate may be enhanced in the minimal 
super-symmetric standard model (MSSM) via intermediate charged Higgs bosons [25 or m 
the Pati-Salam model of quark-lepton unification Similarly, B~^ C^vi may be mediated 
by scalar super-symmetric particles in R-parity violating extensions of the MSSM 0]. In this 
paper, we concentrate on the £ = e and /x modes. The most stringent current upper hmits 
for these modes are B{B+ -> e+z/g) < 1.5 x 10"^ [ij and B{B+ fi+u^) < 2.1 x 10"^ 
at the 90% confidence level. 

It is natural to extend the B~^ — > i'^i'i searches to the corresponding radiative modes. 
Because of the additional photon, helicity suppression does not occur in these modes; hence 
B^ i^i'£'~f decays are predicted to occur with rates comparable to or possibly larger than 
B^ fi^i^n decay. The B~^ —>■ i'^i'e'j decay has been of theoretical interest as a means 
of probing aspects of the strong and weak interactions of a heavy quark system P-jl^. 
In theory, there are two contributions to B^ —>■ i^uej: the internal Bremsstrahlung (IB) 
process and the structure-dependent (SD) process In the IB process, a photon is emitted 
in either the initial or final state. The amplitude for this process is, however, suppressed 
by both helicity conservation and the electromagnetic coupling to the photon. In the SD 
process, the photon is produced in the transition of the spin-0 B~^ meson to a spin-1 off- 
shell vector (B*^) or axial-vector {B') meson. Helicity suppression is avoided in this spin 
exchange. In a recent study by Korchemsky, Pirjol, and Yan (KPY) 12], the predicted 
branching fraction is in the range (2 — 5) x 10~^. Current upper limits on these modes are 
obtained by CLEO: B{B+ e+u^'j) < 2.0 x 10"^ and B{B+ < 5.2 x lO'^ ^ . 



In this paper, we present a preliminary search for the leptonic and radiative leptonic B 
meson decays B^ — > /U+z/^, B^ e^^'eT, and B^ C-'^i'n- A brief description of the Belle 
detector is given in the next section. The analysis procedure and results are described in 
Section ITTTl for B^ — > i^^f^ decays and in Section |TV| for B^ i'^i^e'j decays. Throughout 
this paper, charge conjugation is implied unless stated otherwise. 
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II. THE BELLE DETECTOR 



The event sample that we analyze for this study corresponds to an integrated luminosity 
of 140 fb"^ accumulated at the T(45') resonance and recorded by the Belle detector at the 
KEKB e~^e~ collider This sample contains 152 million BB pairs. We also use an 

event sample of 15 fb~^ integrated luminosity recorded 60 MeV below T(4S') resonance for 
background determination. 

The Belle detector is a large-solid- angle spectrometer based on a 1.5 T superconduct- 
ing solenoid magnet. Tracking and momentum measurements of charged particles are done 
with a 3-layer double-sided silicon vertex detector (SVD) and a central drift chamber (CDC). 
Identification of charged hadrons is provided by a combination of three measurements: spe- 
cific ionization (dE/dx) in the CDC, photon yield in the aerogel threshold Cerenkov counters 
(ACC), and time-of- flight information from a cylindrical array of 128 scintillation counters 
(TOF). Photons are detected in an electromagnetic calorimeter (ECL) system made of an ar- 
ray of 8736 CsI(Tl) crystals surrounding the TOF system. Electrons are identified based on 
the dE/dx measurement in the CDC, the response of the ACC, and the position, shape and 
energy of the electromagnetic shower registered in the ECL. Electron identification efficiency 
is over 95% in the momentum range of this analysis; the fake rate is below 0.5%. Muons are 
identified by a resistive plate chamber system (KLM) located within the solenoid's external 
return yoke. Muon identification efficiency is approximately 90% in the momentum range 
of this analysis; the fake rate is approximately 1%. The Belle detector is described in detail 
elsewhere ll 711. 



III. SEARCH FOR B+ ii+v^ 

The main signature for the B^ f^~^^fj. decay is a single muon having approximately 
2.6 GeV/c momentum. The muon candidate is selected from the charged tracks that are 
produced near the interaction point (IP), and identified by its penetration depth and profile 
in KLM. We use only well-identified muons detected in the barrel region of the Belle detector, 
where muon identification efficiency and fake rate are well understood. Muon candidates 
that are consistent with being produced via a J/ijj ^ fi~^fi~ decay are vetoed. 

To exploit the large missing energy and momentum of the undetected neutrino in the 
B^ ^'^I'fj, decay, we require that there should be only one identified lepton in an event, 
and that the direction of missing momentum be within the detector's fiducial volume 
(|cos6'miss| < 0.88, where O^i^s is the polar angle of the missing momentum vector rela- 
tive to the z axis that is opposite the positron beam line). These requirements make it likely 
that the signal neutrino is the only undetected particle in the event, so that the missing 
mass is consistent with zero. We require —4.0 < M^j^^ < 4.3 (GeV/c^)^. 

The detected particles in the event other than the signal muon candidate are required to 
be consistent with a B decay hypothesis. This is tested using two kinematic variables, Mbc = 
\/-^beam " (Ei PiT and AE = J2iEi- -Ebeam whcrc -Ebeam IS the beam energy measured in 
the center-of-mass (CM) frame and the summations are over all detected particles except the 
signal muon. We require each event to lie within the region 5.1 GeV/c^ < Mbc < 5.3 GeV/c^ 
and —2.0 GeV < AE < 1.5 GeV. The formation of this second B meson's four-momentum 
permits us to determine the reference frame in which the parent of the signal muon is at 
rest ("5 rest frame"). We reject events for which the muon momentum in the B rest frame 
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FIG. 1: Distributions of the event shape Fisher discriminant for signal and background MC samples. 
We require T > 0.1. 

lies outside the range 2.58 GeV/c < < 2.8 GeV/c. 

The dominant background arise from B X^iug decays and from the continuum process 
e~^e~ qq {q = u, d, s, c). There is some additional contamination due to misidentification 
of a hadron as the signal muon. 

The lepton from B — > X^lvi decays (but not the more abundant B — > Xclvi decays) can 
have sufficient energy to populate the acceptance window. These events are suppressed 
by the M\^c^AE cuts, since the second B meson inherits the additional particles associated 
with the Xu hadronic system. 

Since a B meson pair produced at the T(45') resonance has very little kinetic energy in 
the CM frame (7/? ~ 0.06), the event topology of BB events is spherical while that of the 
continuum events, where the light quarks are produced with large kinetic energy, tends to be 
two-jet like. To suppress continuum background, we ex plo it these event shape differences by 
combining several Fox- Wolfram event shape moments fisj into a Fisher discriminant JF 
whose coefficients are trained using Monte Carlo (MC) simulated event samples to maximize 
the distinction between signal and background. (See Fig. [TJ) We also make use of the angle 
Ot between the muon momentum and the thrust axis of the second B meson in the CM 
frame. We suppress continuum background by requiring JF > 0.1 and Icos^tI < 0.56. 

Figure El shows the distributions for signal and background MC events that satisfy all of 
the above criteria except the p^ cut. Signal yield and background contamination are assessed 
in the two-dimensional subregion Mbc > 5.26 GeV/c^ and -0.8 GeV < AE < 0.4 GeV. The 
definition of the signal region, as well as the requirements on p^, event-shape cuts and M^j^g 
are optimized using signal and background MC samples. Using a signal MC sample, the 
signal detection efficiency is determined to be (2.9 ± 0.1)%. 

To estimate the background in the signal region, we define a side-band region by 5.1 < 
Mbc < 5.24 GeV/c^ and —2.0 < AE < 1.5 GeV. This and the signal region are indicated 
in Fig. 121 We fit the Mbc projection of the side-band distribution of the MC background 
sample to an empirical threshold function ("Argus function") J20|. Then we estimate the 
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FIG. 2: Distributions of the muon momentum in the B rest frame satisfying all selection criteria 
except the cut on p^. We require 2.58GeV/c < pj^ < 2.8GeV/c. 
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FIG. 3: Distributions in Mbc-AE^ plane for (a) signal MC events and (b) background MC events 
satisfying all other event selection cuts. The signal region and side-band region are indicated with 
rectangles inside the plot. 

background in the signal region by fitting the Mbc projection of the data side-band with the 
same Argus function where the fit parameters except for the overal normalization are fixed 
to those obtained in the MC fit. The expected background in the signal region, obtained 
by extrapolating the Argus function, is A^bkg = ^'^■'^-I'.t events. Assuming no signal events 
and based on a Poisson fluctuation of background events, we calculate mean expected upper 
limit as 1.7 x 10~^, where backgrounds are subtracted. 

Figure m shows the Mbc distributions for data, background MC, and signal MC events after 
applying all selection criteria except the Mbc cut. We find 4 events in the Mbc signal region. 
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FIG. 4: The Mbc distributions for data (points with error bars) as well as signal and background 
MC events (histograms) for events satisfying all selection criteria except the cut on M^c. The 
signal region is Mbc > 5.26 GeV/c^. 

Without subtracting background, we convert this observation and the above expectation for 
iVbkg into an upper Hmit of 7.99 on the signal yield at the 90% confidence level. 

The systematic uncertainty in the — > analysis includes uncertainties in the signal 
detection efficiency and in the background estimation. Since we do not subtract background 
in this analysis, we consider only the uncertainty in the efficiency. 

The uncertainty in is 0.3%. The uncertainty in track- finding efficiency is 2% per 
track, and that in muon ID efficiency is 2%. The effect of other event selection cuts is 
studied by analyzing a calibration sample of fully reconstructed — > i)(*)°7r"'" decays, 
where we treat the pion as a signal muon and the D^*^° as the accompanying neutrino. 
We compare the efficiency of the event selection procedure between data and MC to infer 
systematic uncertainty in the event selection procedure. The fractional difference of efficiency 
between data and MC is 5.8%. Including the systematic uncertainty, the signal efficiency is 
(2.9 ±0.2)%. 

For a conservative upper limit, we reduce the signal efficiency by its error and obtain the 
following upper limit for the f^^i^fi branching fraction at 90% CL: 

B{B+ /i+z/^) < 2.0 X 10^6 
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IV. SEARCH FOR B+ 



We search also for the three-body radiative decays —>■ e^Vel and —>■ The 
presence of the photon invahdates the monochromaticity constraint on the lepton momentum 
in the B rest frame. Nevertheless, theoretical calculations!^ predict the lepton momentum 
spectrum to peak at the high end of its kinematical range. While the photon's energy is 
typically much softer than the lepton's, it exceeds that of most photons produced at T{4S) 
resonance energy 

We select the primary muon or electron candidate with CM momentum between 1.8 
and 2.8GeV/c using the particle identification information described in section |nj and the 
primary photon candidate as the most energetic of the neutral clusters in the ECL (i.e., 
those not matched to any charged track) above 1.0 GeV that does not appear to be affiliated 
with a 7r° meson. (A photon is considered to be a 7r° meson daughter if, when combined with 
any other photon of energy 100 MeV or more, it gives an invariant mass within ±15MeV/c^ 
of the nominal 7r° mass.) As before, we require that there be only one identified lepton per 
event, that the direction of the missing momentum be within the detector's fiducial volume, 
and that the squared missing mass fall in the range —2.0 (GeV/c^)^ < M^j^g < 0.8 (GeV/c^)^. 
In addition, we require the charge sum of all detected particles be or ±1. 

As before, we calculate Mbc and AE for the remaining particles in the event; here, can- 
didates must satisfy S.llGeV/c^ < Mbc < 5.29GeV/c2 and -0.5 GeV < AE < 0.125 GeV. 
We also find that the lepton momentum (p^) spectrum are different for different sources of 
backgrounds: the lepton spectrum of BB MC events is peaking near 2.0 GeV/c and drops 
out rapidly beyond ~ 2.3 GeV/c while that of the continuum MC sample is much flatter 
than the others and reaches to higher momentum regions (~ 3 GeV/c or beyond). The 
signal yield is obtained by fltting the two-dimensional p£-Mbc distributions. 

Continuum background is suppressed by a cut on the normalized second Fox- Wolfram 
moment (rather than on the Fisher discriminant) of R2 < 0.35 for e^u^'j {R2 < 0.32 for 
fi'^Ufj^'y). The relative fraction of each background component {BB, Xuii'e and continuum) 
are determined by fltting the lepton momentum distribution within the Mbc side-band region 
of 5.11 GeV/c^ < Mbc < 5.20 GeV/c^ to the sum of the three empirically determined shapes 
for these backgrounds. 

For the B^ f^^'^fil mode. Fig. El shows the distributions within the signal region of 
the squared missing mass and photon energy for data and combined background MC (with 
relative weights equal to those measured in the Mbc side-band region). In both cases, the 
shapes of the data and combined background are in good agreement. Similar agreement is 
seen for the B^ e^Vel niode. 

Before fltting for the signal yield, the procedure is tested on an MC sample made from 
a mixture of BB, continuum and B — > Xuii^e processes but without including any signal 
events. The signal yield obtained are consistent with zero for both modes, and we take 
the corresponding upper limits calculated without systematic uncertainty, 1.4 x 10~^ for 
B^ e^Vel niode and 1.0 x 10~^ for B^ as the search sensitivity for the 

corresponding mode. 

We perform a binned maximum-likelihood flt to data in the p^-Mbc plane; the two 
parameters in this flt are the number of signal and combined-background events. Fig- 
ure 10 (Fig. [7|) shows the projections of the flt onto the M^j^g and axes for the 
B^ e^Pel {B^ — > /x+z/^7) mode. The projection onto Mbc is made after applying a 
cut 2.2GeV/c < pi < 2.8GeV/c, and the projection onto pi is made after applying a cut 
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FIG. 5: Data and combined MC distributions of (a) squared missing mass and (b) photon energy 
for the //+f^7 mode after apphcation of ah event selection criteria except the cut on the 

displayed quantity. 

5.26 GeV/c^ < Mbc < 5.29 GeV/c^. 

Since neither mode shows a significant excess above the expected background, we set 
upper limits by integrating the likelihood as a function of the signal yield. The dotted 
histograms in Figs. IHl and [7| show the 90% (statistical error only) confidence level (CL) 
upper limits, added to the background. The fitted yields after efficiency correction are 
970I55Q events in the e^i^el niode, and SOtsgJ] events in the y^^^nl mode. 

Using the statistical error only, the corresponding 90% CL upper limits on the branching 
fractions are 1.2 x 10~^ for 5+ e^Vel, and 0.7 x 10~^ for B^ l^'^'^t^l- 

Systematic uncertainties related to the signal efficiency are estimated as in the B~^ — >■ 
/i"^z/^ analysis. To be conservative, we reduce the signal efficiency by the systematic error in 
efficiency. 

The systematic uncertainty in the yield extraction has three categories: fitting method, 
background estimation, and signal modelling. The uncertainty due to the fitting method is 
estimated by repeating the fit under various scenarios. The systematic effect of background 
uncertainty is estimated by repeating the fit with each background component fraction 
changed by one standard deviation. Similarly, the uncertainty due to signal modelling is 
studied by varying the signal decay parameters in the signal MC generation which is based 
on the KPY model |il2|]. In each category, we take the largest deviation from the default 
procedure as an estimate of systematic uncertainty. Then the three errors are added in 
quadrature for systematic uncertainty in the yield extraction. 

Assuming a Gaussian distribution for the systematic error, we combine the systematic 
uncertainty in the yield extraction to the statistical error by smearing the original likelihood 
function obtained in the signal fitting with a Gaussian function based on this error. Then we 
use this convolved likelihood function to calculate the upper limits. The results, including 
the reduction of efficiency by the corresponding systematic uncertainty, are: 

B{B+ e+z/e7) < 2.2 x 10"^ 
B{B^ ^i-^u^-i) < 2.3 X 10"^ . 
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V. SUMMARY 



In summary, we have searched for — > /^"'"t'ju and — > i'^ug'j decays and found 
no significant excess in any mode. Therefore, we set prehminary upper hmits at the 90% 
confidence level on the branching fractions of B{B^ fJ-^J^iJ.) < 2.0 x 10~^, B{B~^ — > 
e+z/e7) < 2.2 X 10"^ and B{B+ l^'^'^fil) < 2.3 x 10"^. The result on 5+ is more 

stringent than the existing limit by a factor of three. The e^z/e7 mode upper limit is an 
improvement from the existing limit by an order of magnitude and the limit on n'^i'^'j is 
also improved by more than a factor of two. 
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FIG. 6: Distributions of (a) beam constrained B meson mass and (b) electron momentum for 
events within the signal region for the — > e"'"^'e7 mode. The projection onto Mbc is made after 
applying a cut 2.2GeV/c < pe < 2.8GeV/c, and the projection onto pi is made after applying a 
cut 5.26GeV/c^ < Mbc < 5.29GeV/c^. The points represent the data; the curves represent the 
projections of the fit components. Specifically, the dashed curves show the background only; the 
solid curves show the background and signal; the dotted curves show an alternate fit with the signal 
component infiated to its 90% CL upper limit. 
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FIG. 7: Distributions of (a) beam constrained B meson mass and (b) electron momentum for 
events within the signal region for the — > ^'^v^^y mode. The projection onto Mbc is made after 
applying a cut 2.2GeV/c < < 2.8GeV/c, and the projection onto p£ is made after applying a 
cut 5.26GeV/c^ < Mbc < 5.29GeV/c^. The points represent the data; the curves represent the 
projections of the fit components. Specifically, the dashed curves show the background only; the 
solid curves show the background and signal; the dotted curves show an alternate fit with the signal 
component inflated to its 90% CL upper limit. 
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